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I. INTRODUCTION 



In the expression for the polarised-neutron beta-decay rate, one of the terms which can be 
written is a triple scalar product of the neutron spin polarization and the momenta of the 
electron and neutrino; its coefficient is referred to as the U D coefficient" PJ. Such a term is 
time-reversal-odd and, in the absence of final-state interactions, a non-zero coefficient would 
indicate the presence of time-reversal violating interactions in the beta-decay vertex. 

An estimate of the D coefficient from CP-violation in the standard model turns out 
to be extremely tiny, D^m < 10 -12 Q, whereas a prediction from, e.g., the minimal 
supersymmetric standard model (MSSM) is reported as -Dmssm ~ 10~ 7 ^ and more 
exotic mechanisms such as leptoquarks could give even larger results Experimen- 
tal measurements of the D coefficient may soon probe these values; recent papers give 
bounds of D — [—0.6 ± 1.2(stat) ± 0.5(syst)] x 10~ 3 from the EMIT experiment 0] and 
D = [— 2.8±6.4(stat) ±3.0(syst)] x 10~~ 4 from the Trine collaboration (6|; for a review of the 
future experimental prospects see, e.g., Ref. 0. Increasingly precise measurements could 
open a new era for studies of physics beyond the standard model. A complication is that 
there are small but non-zero corrections from the electromagnetic final-state interaction [l[ . 
The correction at the one-photon loop order vanishes in the zero-recoil approximation, and 
(relaxing that approximation) has been estimated by Callan and Treiman as Dps — 10~ 5 0. 1 
If new physics which contributes at or below this magnitude is to be explored, it would be 
desirable to know this value to an accuracy of 1% or less. As previously established in many 
other processes, the ideal tool for such work is low-energy effective field theory (EFT) flo| . 



In this work, we calculate the D coefficient employing a pionless heavy-baryon EFT based on 
heavy-baryon chiral perturbation theory (HB^PT) with photon fields [ll, 12|, T3 |. HB^PT is 



a low-energy effective field theory of QCD and has a systematic expansion scheme (counting 
rules) in terms of small external momenta, symmetry-breaking terms (proportional to the 
pion mass, m n ), and the number of loops. "Small" is in relation to the scale A x rs m p , 
which characterises the physics which is not explicitly included. "Heavy-baryon" indicates 
the approach in which the nucleon mass M N is also treated as a "large" scale, M N m 
A x . Hence the usual expansion parameter of HBxPT is roughly of order m n /M N . (For a 
recent review see Ref. |l4|].) However, because the typical energy of neutron beta decay is 
so small compared even to m n , the pion itself can be integrated out and is not included 
as an explicit degree of freedom. The counting rules for neutron decay in the pionless 
EFT have been worked out by Ando et al [TBI]. We employ a/2n (which governs radiative 
corrections associated with photon loops), and Q/M N (from the heavy-baryon expansion) 
as our expansion parameters, where a is the fine structure constant and Q is a typical 
momentum of neutron beta decay, Q ~ m n — m p — m e (m n , m p and m e are the neutron, 
proton, and electron mass, respectively). These two expansion parameters are numerically 
almost the same, a/2ir ~ Q/M N ~ 10 -3 . 

As we will show the Callan- Treiman result is reproduced in the EFT framework as the 
(D(aQ/M N ) contribution, 2 but the counting suggests that 0(a 2 /2ii) contributions could be 
as large. Such contributions coming from repeated Coulomb final-state interactions have 



For a historical review, see Ref. Q. This result has been confirmed by Refs 16, 17, 3] 
Since, as we will sec, the D coefficient requires a loop integral with an imaginary part, contributions to it 
will be generically a factor of 2tt larger than other terms at the same loop or 1/M N order, which we will 
keep track of: hence 0(aQ/M N ) rather than (D(aQ /2irM N ). 
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already been shown to vanish for zero recoil lq . but radiative effects have not. However 
we are able to show that in the limit M N — > oo all electromagnetic contributions to the D 
coefficient vanish, and so corrections to the Callan-Treiman result must be higher order. 

In fact the counting rules discussed above, while valid for the first three orders, must be 
modified at the order at which the counterterms encode information about the pions which 
have been integrated out. In fact the leading higher-order correction to the D coefficient, 
which we calculate, comes from the pseudoscalar form factor of the weak nucleon current, 
which is enhanced because its scale is governed not by the nucleon mass but by the pion 
mass This allows us to conclude that the residual error is of the order of 1 part in 10 3 , 
and certainly smaller than 1%. 



II. THE D COEFFICIENT AND EFFECTIVE LAGRANGIAN 

The general expression for the differential neutron decay rate dT is well known for the case 
wherein only the neutron is polarized: 

dT ~ (G F V ud ) 2 3q 2 W i e E 2 

dE e dn p dQp„ - (2tt)« ^+^9 A )\Pe\^eK 

1 , Pe ' Pv . a ( A Pe . u Pv . n Pe X p u \\ 

1 + a — hn- A— + B— + D , (1) 

E e E u \ E e E v E e E v J J 

where Gf is the Fermi constant, V u d is a CKM matrix element, and g A is the axial-current 
coupling constant. Here E e and p e (E u and p v ) are the electron (neutrino) energy and 
momentum, n is the neutron spin polarization vector, and a, A, B, D are the correlation 
coefficients. The standard lowest order expressions for the correlation coefficients are 

1 + Zg\ 1 + Zg\ 1 + 3gl y 1 

We can reproduce these expressions for the coefficients in the leading order (LO) EFT 
calculation. In addition, at NLO the 0(a/2n) radiative corrections and the 0(Q/M N ) 
recoil corrections (including the weak magnetism term) to the decay rate T and correlation 



coefficients a, A, and B have already been reported in Ref. [15]. A non-zero contribution to 
the D coefficient, due to the electromagnetic final-state interaction, has been obtained by 
Callan and Treiman from one-photon-loop contributions and we will see that the same 
expression comes out at NNLO. 

The standard chiral counting rules for HBxPT with and without photons are discussed in 



Refs. [11|, [12|, [13j. However as mentioned above, and discussed in more detail in Ref. [15 
the smallness of the energy released in neutron beta decay modifies the counting so that 
we have two parameters, Q/M N ~ Q/A x and a/2n, which are is numerically of the same 
magnitude. Thus we treat them as a single common expansion parameter in defining LO, 
NLO, etc. 

The effective Lagrangian for the calculation of the D coefficient in neutron beta decay 
reads 15 



£/3 — £ewy + £7V7V7 + £-NNeu , (3) 

where C eui is the standard Lagrangian for electron, neutrino, and photon, Cnn^ is the heavy- 
baryon Lagrangian for nucleon interaction with a photon up to NLO (the 1/M N order), and 
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£-NNev is the V-A interaction Lagrangian between the nucleon and lepton currents up to 
NLO: 

£ e „ 7 = -- A F af3 F a p - -L(«9 . A) 2 + &(i 7 • D - m e )^j e + $„i 7 • , (4) 
c NNl = NUvDN 

+ 2AT N N '( {V ' D? -° 2 - 2^°' + 5A*sTr[/+,]))jV, (5) 

C NNev = -^^4 7q (1 - 75 )^(iVV + b Q " Zg A S a ]N 

+ 2i^~ iVtr+ " 9 " P ~ 2gAvasP Kd ~ ^ ~ 2t ^[S a , S-(d + a)]) iv) , 

(6) 

where = d a Ap — dpA a and _D Q is the covariant derivative of QED. £a is the gauge param- 
eter and we choose the Feynman gauge £4 = 1. v a is the velocity vector with the condition 
v 2 = 1. In the rest frame of the neutron, = (1,0) and 2S* M = (0, a). Furthermore, the 
photon couples to the nucleon via the charge operator Q = |(1 + r 3 ), giving = 2QF Q/ g 

and f + = f + — |Tr[/ + ]. /i v and fi s are isovector and isoscalar nucleon magnetic moments; 
/V = /Up - yUn = 4.706. 

The counting rules discussed above implicitly assume that higher-order terms in the La- 
grangian are suppressed by successively more powers of A x w M N . In reality however the 
break-down scale of the EFT is governed by the lightest degree of freedom which has been 
(at least conceptually) integrated out. Here, this is the pion, which is substantially lighter 
than A x . The low-energy constants (LECs) of the pionless theory may therefore be governed 
by inverse powers of the pion mass, and hence larger than expected. 

To the order we are working, the only LECs which enter are coupling constants such as g A 
and the magnetic moments, and these are fixed by their well-known experimental values. 
At higher orders though, potentially l/m^-enhanced form-factor terms start to show up. In 
fact at the next order there is only one such term, the induced pseudoscalar form factor of 
the nucleon. This is discussed in the final section. 



III. AMPLITUDES 



The relevant Feynman diagrams for the NNLO calculation are shown in Figs. [H [2], and[3j A 
contribution to the D coefficient is obtained from the interference between the imaginarypart 
of an amplitude from a loop diagram and the real amplitude of a tree-level diagram [if, 0] . 
Fig. [1] contains the tree diagrams and those loop diagrams which have an imaginary part and 
which contribute to the D coefficient; these together give the 0(aQ/M N ) contribution. The 
others do not contribute, in the case of Fig. [2] because they have no imaginary part. The 
diagrams of Fig. [3] (which with the leading tree graph would give 0(a 2 /2ir) contributions) 
will be considered in section 5. There is no 0((Q /M N ) 2 ^ contribution since this comes from 
tree amplitudes only. 

We first consider the tree-level amplitudes. The LO amplitude from the diagram (a) in 
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FIG. 1: Diagrams for the 0{aQ /M N ) contributions to the D coefficient. In all diagrams the LO 
nucleon vertices are denoted with a dot, and NLO ones (0(Q/M N )) with a cross. A wavy line 
denotes a photon. Diagrams (a) and (b) are tree-level LO and NLO diagrams. Diagram (c) is NLO 
(0(q27t)) while diagrams (d), (e), and (f) are NNLO as they include 0(Q/M N ) corrections. 




(k) (1) (m) 



FIG. 2: Diagrams for the 0(aQ/M N ) amplitudes which do not contribute to the D coefficient. 
See the caption of Fig. Q] for the details. 



Fig. [T] reads 



G V 

M {a) = -^^u e (p e )l a {l - l->)v v (p u )xliv a -2g A S a )xn 

GFVud M Q , (7) 



V2 

The 0(Q/M N ) tree-level amplitude from the diagram (b) in Fig. [1] reads 

x^xj {(9 aP - v a v? + 2g A v a S?)(p p + p n ) p + 2MS°, S ■ ( Pp - p n )\) X n , (8) 

where p° and p" are the residual four-momenta for proton and neutron, respectively. We 
define them as p° p = P" p — m p v a where ?° and P p are four-momenta for neutron and 
proton, respectively: P% = m\ and P p = m p . 




FIG. 3: The non-trivial diagrams for the 0(a 2 /2n) amplitudes. See the caption of Fig. [T] for 
details. 



We now consider the loop diagrams in Fig. [TJ The diagram (c) in Fig. [TJ has already been 
calculated in Ref. 15 



giving 

M (c) 



^ GFVud (M F + M 1 (-m e f 2 )) 



V2 



with 



Mi = M e (p e )^7 Q (l - 75)u„(p„)xJ(ua - 2g A S a )Xn , 
F = f 1 + 2E e (f 2 -f Q ), 

where fo, fi, and f 2 are loop functions whose definitions are given in Appendix A. 
From the diagrams (d), (e), and (f) in Fig. (TJ we have 



M t 



(d) 



GfV u <i e 2 



We7 M (^ e -/ + "ie)7 a (! - 75)^ 



Mi 



(e) 



y/2 2M N J l vll 2 (l 2 -2p e -l) 

x xl ((flV ~ V")( f + 2p p ) v - 2/i p [S^, S ■ I}) (v a - 2g A S a )xr, 
G F Vud e 2 fu e ^(^> e -/ + m e )7 a (l - j 5 )v v 



^2 2M N J l vll 2 (l 2 -2p e -l) 
x xl(g a p - v a v p + 2g A v a Sp + 2fi v [S a , Sp\)(l + p p + p n ) P Xr 



Mi 



G F V ud e 2 



if) 



2M W J Z (vl) 2 l 2 (l 2 -2p e -l) 



Z 2 - / 2 - 2p p • /) 



with 



(2vr 



|4 ' 



(9) 



(10) 
(11) 



(12) 



(13) 



(14) 



where we have used the fact that the nucleon kinetic energy is 0(Q 2 /M N ) to set v ■ p v ~ 
in the nucleon propagator. 



IV. THE D COEFFICIENT FROM HEAVY-BARYON EFT 



In calculating the decay rate, which is proportional to the modulus of the amplitude squared, 
we are taking the product of one diagram and the time-reversed version of another, e.g. 
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M( a )M* d y Considering the same pair in the opposite order, M^MT a y it is obvious that the 
time-reversal odd term in ■ (p e x p v ) will have the opposite sign and a complex-conjugated 
coefficient compared with the original order. Thus when we sum the two, we can only get 
a non-zero result if the coefficient has an imaginary part, and this must come from a loop 
integral with a physical intermediate state. 

From the interference between the LO amplitude from the diagram (a) and the imaginary 
part of the leading loop diagram (c), we have as the potential NLO contribution to the D 
coefficient 

D 

D<x^2 (M (c) M ( * a) + M (0) M ( * C) ) = , (15) 

spin 

where the superscript D on the summation denotes that we sum over the spins of the the 
electron, neutrino, and proton but not that of the neutron, and keep only terms contributing 
to the D coefficient, those proportional to n-{p e xp u ) where n = Xn&Xn- (The corresponding 
4- vector spin polarisation N* 1 is given by = xli^S^Xn, and ie^ vafi ] v ^N v p ea p v p = n ■ (p e x 
p u )-) This NLO term vanishes because there is no such term in ^ spin -MoM 0) i, an d not 
through details of the loop integrals Fq and f\. Thus the D coefficient has no contribution 
up to NLO. 

We now consider 0(aQ/M N ) interference terms. From the amplitudes from the NLO 
(0(Q/M N )) diagram (b) and leading loop diagram (c), we have 

D 2 

(M (C) M ( * 6) + M (6) M ( * C) ) = -4(G F V ud ) 2 ^-e^v,N uPeaPu p 

spin 

x {2ImF (E e (-g A n v + g A + g 2 A - /v) + K (-g A Hv + 9a~ g A + fJiv)) 

+2lmf 2 m 2 e (g A fi v - g A - g\ + fi v )} . (16) 

From the LO amplitude Mr a ) and the NNLO (0(aQ/M N ) amplitudes M^^j), we have 

D 2 

(M id) M( a) + M (a) M ( * d) ) = -4(G F V ud ) 2 ^-e^%N u p ea p u p 

spin 

x(l - g A ) (-%g A E v lmh + 2^(1 + 3g A ) (2Im/ 4 + E e lmf 5 + m e 2 Im/ 6 )) , (17) 

D 2 

(M {e )M{ a) + M (a) M ( * e) ) = -A{G F V ud ) 2 ^—e^ Vll N v p ea p v p 



2M N 

spin 



{ x2ImF (E e (g A /j, v - g A - g 2 A + fi v ) + E u {g A ^ v - g A + g 2 A - /i v )) 
+2(g A fi v -g A -g 2 A + /i v ) {-E e lmf 5 + (-m 2 e + 2E 2 )(lmf 2 - Im/ 6 )) 
+4(/i v - g A )lmU} , (18) 



D 2 



J2 {M (J) M* a) + M (a) M* u) ) = -A(G F V ud ) 2 —e^v^N u p ea p vP (-8^(1 - g A )E„Img 4 ) , 

spin 

(19) 

where explicit expressions for the loop functions and g± are given in the Appendix. 
(Throughout, we have used unit normalization for the heavy-baryon nucleon spinors). 



s 



To obtain the D coefficient, we have to divide out the common factor 4(G fV uc i) 2 E e E u (l+3g\) 
which also appears in the angle-independent term in the leading matrix element squared, 

Espin^a)^)- 

We thus have the final expression for the 0(aQ/M N ) contribution to the D coefficient for 
neutron beta decay: 

Do? = rr^2^^(( 1 + 3 ^)((^-^)- 3 Mi-^)) 

m 2 \ 
+ ((3 + g A ){fi v - g A ) + 3/i p (l - g A )(l + 3g A ))j . (20) 

where (3 = \p e \/E e . This is exactly the result obtained by Callan and Treiman [§[. 



V. TERMS HIGHER-ORDER IN a 

Since, as discussed above, the two scales a/2ir and Q/M N are of the same size, 0(a 2 /27r) 
terms, if they exist, could be as large as the well-known 0(aQ/M N ) contribution which we 
reproduced in the previous section. We therefore need to check these too. The relevant 
diagrams are shown in Fig. [3], and their amplitudes are as follows: 



M {A ) oc 



Ue{Pe)j> (A + jj + m e) j> ty e +/ + m e ) jjj. - 7sK(P»/) f , > 
Vlvkk 2 (k-l) 2 (l 2 -2Ve-D(k 2 -2p P -k) ^ ^ 



M(c) oc 



J k vlvkk 2 (k-l) 2 (l 2 - 2p e ■ I) (k 2 - 2p e ■ k) 
M(B) * JJ k vlv.(l-k)k 2 (k-l) 2 (l 2 -2p e .l)(k 2 -2p e .k) X > ~ 29aS )Xn 

U e (Peh a {i> e + ¥' + ™> e ) j) {jf> +/ + ^ + m e ) 7 CT +/ + 77l e ) 7 a (l - 1s)v v {p v ) 



I Jk 



vll 2 k 2 (I 2 - 2p e ■ I) (k 2 - 2p e ■ k) {{I + k) 2 - 2p e ■ (I + k)) 



xxUv a -2g A S a )x, 



M {D) oc / / ^^'^^L"^^ ^ - 2*5")*. • (21) 



J k vlv(l + k)vkl 2 k 2 (l 2 -2p e -l) 

Clearly the structure of the amplitude M^d) is the same as the leading one-photon-loop 
diagram Mm, though the integral is different. So interference with M cannot give a con- 
tribution to the D coefficient. The other three cases are more complicated (though it is 
worth noticing that the tensor structures of Mmi and M^b) are the same; again only the 
integrals differ). However explicit calculation — without the need actually to calculate any 
scalar two- loop integrals — again shows a vanishing contribution to the D coefficient. (The 
calculations are most efficiently done using e.g. the package "Tracer" on Mathematica (20| : 
we do not give details of the reduction to scalar integrals as in the Appendix because they 
are not required.) 

Having shown through explicit calculation that there are no two-photon-loop contributions 
to the D coefficient in the heavy-baryon limit, it is interesting to consider if this can be 
generalized. The crucial features which make it possible are that in this limit, as can 
be seen in the expressions for the diagrams (c) and (A-D), the amplitude factorizes into 
a hadronic and leptonic part. Since the leading-order photon coupling to the proton is 
spin-independent, no matter how complicated the photon exchanges and dressings are, the 
hadronic tensor is unchanged. Furthermore, since the heavy-baryon propagator does not 



9 



depend on the three-momentum of the proton, the loop integrals can only pick up factors 
of the electron momentum and the velocity vector, limiting the complexity of the structures 
which can appear in the leptonic tensor. Thus it is plausible that many-photon effects do 
not in fact generate new structures. Recall that the vanishing of the D coefficient in the one- 
and two-photon-loop case is due to the tensor structure, and not to details of the integrals 
involved. 

In fact it was shown long ago that one class of corrections vanish to all order in a, namely the 
repeated exchange of Coulomb photons [l6|. In the heavy-baryon limit these only introduce a 
phase shift in the final-state wavefunction, and cannot induce a D coefficient. However in the 
relativistic theory crossed photons, or multiple overlapping dressings of either fermion line, 
cannot be accounted for so easily as the interaction kernel becomes arbitrarily complicated. 
In the heavy-baryon formalism however the complexity is much reduced, as indicated above. 
As we will demonstrate in what follows, we can show that no corrections survive at any 
order. 

The amplitude resulting from the exchange of any number of (non-magnetic) photons 
(crossed or not) between the electron and proton, and from photon-loop dressing of one 
or more photon-electron or photon-proton vertices, has the same general form as the single 
photon diagram (c), albeit with multiple integrals and many insertions along the electron 
line. The basic structure can be written schematically as follows: 

M X OC ( [ . . . [ *e(p e )6({fr})7 a (l - 7^M_ ) U y a _ 2 Sa) (22) 

\J h J ln vl>...lf...(l> k 2 -2p e .l> k )...J Xp{ )Xn ' 1 ' 

where the l[ are linear combinations of the loop momenta k, and 0{{l{\) is composed only 
of the building blocks m e , {^}, $ e , f and j a . . . 7 " (the last being for a photon loop dressing 
one or more vertices, and with . . . standing for more of the same). 

When the loop integrals are performed, all loop momenta If in the numerator either become 
pa Qr v a^ Qr a p a j r gj ves ga/3 _ ^hus a ft er integration, the whole structure represented 

by 0({li}) can only give a small number of terms, namely /, ^ and [^ e , f\, and the same 
multiplied by 75. Furthermore since 7s7 Q (l — 75) = 7 a (l — 75), the structures with 75 in 
them are redundant. Hence the bottom line is 

4 

Mx = Yl J « K, v ■ Pe )u e (p e )O nla (l - 75 K(p,)4K - 2g A S a } X n , (23) 

n=l 

where the integrals I* will depend on the particular graph we are considering, and the 
operators O n are the four listed above. 

When we calculate the decay rate from the amplitude we need expressions like M^My, with 
a sum over the spin of the proton, electron and neutrino. Using completeness relations such 
as ^ s u e (p e , s)u e (p e , s) = j) e + m e , we end up with terms like the following: 



e 
4 

Tr[ 7a (l - l5 )O m (k + m e )6 nlp {l - 75)^] x 

m,n=l 



xi(v a -2g A S a )(v^-2g A S^ Xn , (24) 

where O m = 'joO^'jo = Finally, we note that O m {j) e + m e )O n is itself just a combina- 

tion of the O m , giving just four structures to be calculated: 

Tr[ 7a (l - 75)0^(1 - 75 )AJ xUv a - 2g A S a )(v 3 - 2g A S p ) X n • (25) 
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Only O n = j) e and ji give non-vanishing results, and they do not generate the structure 
e af3cTTPePv v ° X n S 7 Xn- So in the heavy-baryon limit, there are no contributions to the D 
coefficient at any order in a. Corrections to the leading aQ/M N result will be 0(a 2 Q /2ttM n ) 
and 0(aQ 2 / M 2 ) — presumably at least another factor of 10~ 3 down. 3 



VI. HIGHER-ORDER CORRECTIONS AND CONCLUSIONS 



In the preceding sections, we have shown that the Callan-Treiman result for the D coefficient 
in neutron (3 decay is the leading non-vanishing contribution in a heavy-baryon EFT in 
which the expansion parameters, a/27r and Q/M N , are both of the order of 10 -3 . The 
non-zero contribution to the D coefficient appears because of the electromagnetic final- 
state interaction at 0(aQ/M N ), whereas there are no contributions from the 0(a 2 /27r) 
and 0[{Q / M N ) 2 ^j terms. The D coefficient is NNLO with respect to the decay rate, and 
hence very small (~ 10 -5 ). Nonetheless there are hopes that this might be experimentally 
accessible in the not-too-distant future. The interest in such experiments would, of course, 
be to detect deviations which might indicate new physics, and to interpret such a result 
the accuracy of the Callan-Treiman prediction must be known. At first glance our results 
suggest that the next correction would introduce a relative error of order 10" 3 , which should 
be small enough to be unimportant for many years to come; in particular this would be 
small enough to allow detection (if sufficiently sensitive experiments could be carried out) of 
the prediction from the MSSM of a contribution to the D coefficient of the order of 10 -7 

However at the next order (N 3 LO) there will be contributions from the third-order La- 
grangian, and these will not only be the 1/M 2 terms required by Lorentz invariance, but 
will also include new structures whose scale is governed by the lightest degrees of freedom 
which have been integrated out — in this case the pion. Potentially therefore the corrections 
to the 0(aQ/M N ) result are 0(aQ 2 /ml), giving a relative error of perhaps 5%. 

Looking at the diagrams of Fig. 1, such insertions from the third-order Lagrangian could in 
principle replace any of the crosses, though that in Fig. If could only give a proton mass 
shift. However because we need to take the imaginary part of the loop graphs, the nucleon is 
on-shell everywhere and most of the vertex corrections are just form-factor corrections to the 
structures we have already considered — e.g. fi p — > /i p (l + (r 2 i . I )q 2 /6) — which are therefore two 
orders down. 4 Only one potentially-enhanced third-order effect remains, and that is where 
the cross in Fig. lb and Fig. le represents the pseudoscalar form factor of the nucleon, the 
leading contribution to which is given by the exchange of a pion between the nucleon and 
leptons. 

This pseudoscalar vertex is easily calculated from HB^PT to be as follows (where (q 2 — m 2 ) 
in the pion propagator has been replaced with — m 2 . since the momentum transfer is very 



3 The referee of this paper has drawn our attention to the paper by Gross 2l| on two particles of unequal 
mass interacting via relativistic exchange of some boson, and in particular to the effective one-body Dirac 
equation obeyed by the light particle when the heavy particle's mass is taken to infinity. It is possible 
that the all-orders vanishing of electromagnetic contributions discussed here may also be explicable within 
that framework. 

4 If we work in HB^PT with explicit pions, these form factors arise partly from pion loops. The leading 
term which is 0{Qm- K / '(4m f£) simply rcnormalises the isovector magnetic moment and the next term in 
an expansion in the photon momentum transfer gives a term of relative size Q 2 /m^. 
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much less than m n ): 

- iMp S = i = -q a S -0-7 (1 - 75) , (26) 

V 2 "\ 

There are two contributions to the -D coefficient, arising from the interference of the ampli- 
tudes (la) and (le), and of (lb) and (lc) (with the cross in (lb) and (le) now representing 
the new vertex). These give a final result of 

Ar = 1 1 2 2^ 1 (-4(E e + ^)Im/ 2 + 2Im/ 5 + 2E e Im/ 6 ) 



1 + 3g% ml 

J ami 9 2 A 2(E e + E u ) 

1 + 3g 2 A ml (3 E e 



(27) 



With the approximation E e + E„ = m n — m p this has the value —5.88(pf ax /p e ) x 10 8 
compared to the leading result of (0.228(p™ ax /p e ) + 1.083(p e /p™ ax )) x 1CT 5 — a correction of 
between -0.5% and -2.5%. 5 

With the only (l/m n ) -enhanced N 3 LO term explicitly calculated, we can now say with 
confidence that our final result for the D coefficient, 

D = D CT + D w (28) 

where Dqt an d D n are given by Eqs (1201) and (I27j) . should be accurate to better than 1%, 
with expected corrections being (9(10~ 3 ). 



Acknowledgments 

SA would like to thank J Behr and H M Shimizu for discussions, H-W Fearing and M Igarashi 
for communications. JAMcG would like to thank M Birse for discussions and for reading 
the manuscript and T Cohen for discussions. The work of SA and JAMcG is supported by 
STFC grant number PP/F000448/1. The work of TS is supported by the Japan Society for 
the Promotion of Science, Grant-in-Aid for Scientific Research(c) 20540270. 

Appendix: Loop functions 

The loop functions for the D-calculation in HB^PT are defined as 



1 



(v ■ I + irj)(P + ir])(P — 2p e ■ I + irj) 
where 

/x 4 " d f d d l 



/o, (A.29) 



5 We have used the current PDG value of g A — 1.2695 ± 0.0029. 22| The error on g A induces a 0.1% error 
in D at p e = p™ ax . 
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and D-is the space-time dimensions, d = 4 — 2e. Furthermore 



[vlP(P-2p e -l) 

vvu + g^U + W e + p»y)u + pX/e , (A.31) 



^./i+tf/a, (A.30) 



/^•/Z 2 (/ 2 -2p e -Z) 

l (v.iyp(p-2 Pe .i) = v " vVgz + ^ + + + • (A ' 32) 

where we have suppressed the irj terms in the propagators in those expressions. To calculate 
the contribution to the D-term, we need imaginary part of the six loop functions in the 
followings, and thus have 

Im/i = -lmg 4 = — - , lmf 2 



Im/ 4 = -L/?£ e , Im/ B = -i-I J l m / 6 = _J_^_. (A.33) 

167T 167TP 167T jjE e 
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